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Abstract — This paper studies synchronization in coupled 
nonlinear dynamic networks with unknown parameters. Adap- 
tation can be added to one or several elements in the network, 
while preserving the global synchronization conditions derived 
in [22], [26]. This implies that new nodes can be added 
to the network without prior knowledge of the individual 
dynamics, and that nodes in an existing network have the 
ability to recover dynamic information if temporarily lost. 
In addition, when the individual elements feature sufficiently 
rich stable dynamics , as e.g. in the case of oscillators, 
then adaptation actually leads to an exact estimation of the 
unknown parameters. Different kinds of "leaders" are also 
discussed in this context - one type of leader can specify overall 
trajectories for the network, while another can concurrently 
specify dynamic parameters. 

I. INTRODUCTION 

While the study of synchronization of coupled dynamic 
systems has a long history [20], [23], it is still an extremely 
active research topic [2], [4], [8], [9], [11], [18], [19], 
[24], [25]. In two recent papers [22], [26], we proposed 
a new theoretical analysis tool, which we called Partial 
Contraction Theory, which is based on Contraction con- 
cept [13]. We used Partial Contraction Theory to analyze 
the collective behaviors of dynamic networks with arbitrary 
size and general connectivity. Synchronization conditions 
were derived for networks with diffusion-like couplings, 
which essentially require the network to be connected, the 
maximum eigenvalue of the uncoupled Jacobian matrix to 
be upper bounded, and the coupling strengths to exceed a 
threshold. 

In this paper, we extend the results in [22], [26] to cou- 
pled networks with adaptation. We show that synchroniza- 
tion still occurs under very similar conditions if adaptation 
is added to one or more than one elements in a network. 
Estimated parameters will converge to real values if the 
stable system behaviors are sufficiently rich or persistently 
exciting, which is the case when the individual elements 
are oscillators, for instance. This result is conceptually new 
as it differs with most existing adaptation models, in the 
sense that the adaptation here only adds to a small part 
of the network, and the target behavior is synchronization 
rather than an explicit desired trajectory. Imagine for in- 
stance that some elements in a network lose the normal 
values of some parameters. They can recover the losing 
information from the rest of the network, as long as there are 
still elements holding this information. The nodes holding 
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the real parameters can be considered "knowledge-based 
leaders". They differ with the usual "power-based" leaders 
or virtual leaders such as those in e.g. [8], [11] and [22], 
[26], which specify desired trajectories for the network by 
unidirectionally coupling to it. With adaptation, we can also 
add new nodes into a network without knowing the dynamic 
properties of the nodes already in the network. 

Section |ll] briefly reviews Contraction and Partial Con- 
traction Theories, and the main results we derived for 
synchronization. Adaptation for two coupled systems is 
studied in Section|III| and the results are extended to general 
coupled networks in Section ffVl Brief concluding remarks 
are offered in Section IVl 

II. Synchronization in Coupled Networks 

A. Graph Theory Preliminaries 

Let us first introduce some basic Graph Theory con- 
cepts [3], [6], [15] which will be used in the rest of the 
paper. 

A graph G is composed of a set of n nodes and a set 
of t links. If there is a direction of flow associated with 
each link, G is called a directed graph, otherwise it is 
undirected. G is connected if any two nodes inside are 
linked by a path. The components of the adjacency matrix 
A(G) G M nx " are defined as [A]y = 1, for an undirected 
graph if there is a link connecting nodes j and i, and for a 
directed graph if there is a link from node j to node i, and 
as [A]ij — otherwise. The valency matrix D(G) G M. nxn 
is a diagonal matrix with [D]« = Yfj=i a ij- The matrix 
L(G) = D — A is the Laplacian matrix, which is symmetric 
and positive semi-definite if G is undirected. Its second 
minimum eigenvalue A2(L) is the algebraic connectivity, 
which is zero if and only if G is not connected. The first 
eigenvalue is always zero, corresponding to the eigenvector 
[1,1,..., If. 

Assign an arbitrary orientation a to an undirected graph 
G. We get the incidence matrix = D(G T ) G R nXT . For 
each oriented link k which starts from node i and ends at 
node j, [D]j/- = 1 and [D]j& = —1. All the other entries of 
D are equal to 0. Moreover, 

L(G) = D(G CT ) D T (G ff ) 

If the graph is weighted, we have the weighted Laplacian 
matrix 

L/c = D K. D T 

where JC G M TXT is a diagonal matrix with the A; th 
diagonal entry [JC]k = K,j corresponding to the weight 
of the k th link. If G R mxm is a matrix, K is block 
diagonal. Similarly D has block entries I, —I and 0. 



B. Contraction and Partial Contraction Theories 

In two recent papers [22], [26], we studied synchro- 
nization behaviors of coupled dynamic networks. Here we 
briefly introduce the theoretical analysis tools. 

Contraction Theory [13]: Consider a nonlinear system 
x = f (x, t) where we assume f (x, t) is continuously 
differentiable. Consider a virtual displacement Sx between 
two neighboring solution trajectories. We have 

— (Sx Sx) = 2 (5x — — Sx < 2 X ma x Sx Sx 
at ox 

where X max (x, t) is the largest eigenvalue of the sym- 
metric part of the Jacobian matrix di /dx . Hence, if 
Xmax (x, t) is uniformly strictly negative, any infinitesimal 
length ||#x|| converges exponentially to zero. By path 
integration at fixed time, this implies in turn that all the 
solutions converge exponentially to a single trajectory, in- 
dependently of the initial conditions. 

More generally, consider a coordinate transformation 
Sz = &Sx where 0(x, t) is a uniformly invertible square 
matrix. We have 

-(Sz T Sz) = 2Sz T (0 + ©!!)©- 1 Sz 
at ox 

so that exponential convergence of ||5z|| to zero is guaran- 
teed if the generalized Jacobian matrix 

F = (0 + ©|i)© 1 
ox 

is uniformly negative definite. Again, this implies in turn 
that all the solutions of the original system converge expo- 
nentially to a single trajectory, independently of the initial 
conditions. Such a system is called contracting . 

Partial Contraction Theory [22], [26]: Consider a non- 
linear system in the form x = f(x, x, t). Consider the 
virtual, observer-like auxiliary system 

y = f(y,x,*) 

and assume that this system is contracting. By construction, 
y(t) = x(t) is a particular solution of the y-system. Thus, 
if another particular solution of the y-system verifies a 
smooth specific property, then all trajectories of the x- 
system verify this property exponentially. Such a x-system 
is called partially contracting. 

C. Synchronization of Coupled Networks 

In this section, we list two main results in [22], [26], both 
of which will be extended through the rest of the paper. 
Two Coupled Systems: Consider the coupled systems 

xi = h(xi,t) + g(xi,x 2 ,i) 
x 2 = h(x 2 ,t) + g(xi,x 2 ,t) 

where xi, x 2 <E K m are the state vectors, and g could be 
any combination of xi , x 2 and t. 

Theorem 1: If function h in is contracting based 
on a constant 0, xi and x 2 will converge to each other 
exponentially, regardless of the initial conditions. 



The proof is based on the auxiliary system 
y = h(y,t) + g(x x (t),x 2 (t),t) 

which is contracting and has two particular solutions y = 
xi and y = x 2 . 

Coupled Networks: Consider a coupled network con- 
taining n elements 

Xj = f (xf , t) + } j Kjj (xj - Xj) i = l,...,n (2) 

where Nl denotes the set of indices of the active links of 
element i, and the couplings are symmetric positive definite, 
i.e., Kjj = > 0.. 

Theorem 2: All the elements within the coupled net- 
work will reach group agreement exponentially if 

• the network is connected 

• X max (Ji S ) is upper bounded 

• the coupling strengths are strong enough 

where Ji S is the symmetric part of the Jacobian af , 
The proof is based on the auxiliary system 

n 

which is contracting for appropriate choices of the constant 
matrix K > 0. All solutions of thus verify the smooth 
specific property x x = • ■ • = x n exponentially. 

In fact, we can express the synchronization condition 
more specifically as 

A m+ i(LAc) > maxA maa; (Ji S ) uniformly (3) 

i 

where is the weighted Laplacian matrix with each 
weight corresponding to the coupling strength of that link. 

In addition, the coupling forces in a network can be more 
general, such as uni-directional, positive semi-definite, or 
nonlinear. 

Note that the condition l|3} is always true if J; s < 
and the network is connected. For example, assuming 
f (xj , t) = 0, system (|2j represents a schooling or flocking 
model [22] which provides an effective distributed solution 
for group agreement problem [8], [11], [18] or rendezvous 
problem [1], [12]. 

III. Two Coupled Systems with Adaptation 

Consider two coupled systems as in Q, but assume that 
a parameter vector a is unknown to the second system. 
To guarantee state convergence, we generate an estimated 
parameter a through an adaptation mechanism. Specifically, 
the x 2 dynamics is replaced by 

x 2 = h(x 2 ,a,t) + g(xi,x 2 ,i) (4) 
a = PW T (x 2 ,£)x, x = xi — x 2 

with constant symmetric P > and W(x 2 ,i) defined as 

h(x 2 ,a, t) = h(x 2 ,a,t) + W(x 2 ,i)a 



with a = a a . A similar adaptive technique was 
used in [13], [14], but is generalized here in the sense 
that the couplings are bidirectional. The system structure 
is illustrated in Figure ^ 
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Note that although the diffusion couplings are only based on 
variable v, full -state f eedback is needed for adaptation in this case. 
See Appendix IVI-AI for the boundedness proof. The simulation 
result is illustrated in Figure [5] D 



Fig. 1. The structure of two coupled systems with adaptation. 

Theorem 3: In system (0}, x converges to asymptoti- 
cally if X! is bounded and h is contracting. 
Proof: Define the Lyapunov function 



V 
V 



1 ~T ~ 

2 (X X 



E^P^a ) > 



x ( h(xi,a,t) - h(x 2 ,a,t) ) 



x T 



1 dh 

<9x 



(x 2 + x*)dx x < 



The boundedness of xi implies that of x 2 . Assuming all 
the functions are smoothly differentiable, the boundedness 
of V can be concluded since all the states including a 
are bounded. According to Barbalat's lemma [21], V and 
therefore x tends to asymptotically. □ 
Note that a also tends to 0. Furthermore, since x is also 
bounded, we have the asymptotic convergence of x to zero, 
which leads to the convergence of W(x 2 ,i)a to zero. In 
particular [21], if 

r t+T 

3 a > 0, T > 0, Vt > / W T Wdr > al 



then a converges to zero asymptotically, too. 

The boundedness of xi is trivial if g = g(xi,t), 
which is a classical observer structure with xi dynamics 
independent. If g = g(xi,X2,t), we have 

x 1 =h(x 1 ,a,<) + g(x 1 ,x 1 - x,t) = e(xi,x, t) 

where x is bounded. Thus the boundedness of xi is deter- 
mined by the Input-to-State Stability [10] of e. 
Example 3.1: Consider two coupled FitzHugh-Nagumo (FN) 
neurons [5], [16], [17], a famous spiking neuron model, 



i>i = wi(a — «i)(«i — 1) — Wi + 1 + k(v 2 — vi) 
Wi = pv\ — "fWi 



i>2 = V2(& — V2)(V2 
W 2 = /3v 2 — JW2 



U>2 + I + kOUl 



V 2 



where v is the membrane potential and / external stimulation 
current. We have x = [ v ffl ] , 
and 



[a I 7 P] 1 >Q 



h(x,a, t) 



v)(v - 1) 
Pv - 



- w + 1 ■ 



2kv 



which is contracting if the coupling gain k is larger than an explicit 
threshold [22]. The adaptive law is thus 

k = PW T 5c 
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Fig. 2. Simulation result of Example 13.11 The real parameters are a = 
5.32, ft = 3, 7 = 0.1 and I = 20 . The coupling gain fe = 15 . The 
matrix P = diag{0.6, 30, 0.002, 0.4} . All initial conditions are chosen 
arbitrarily. The plots are (a) states versus time, (b) estimator error a 
versus time. 



IV. Adaptation in Coupled Networks 
A. Basic Results 

Consider now the coupled network (0. We show that 
similar synchronization conditions as those in Theorem [2] 
can be derived if adaptation is added. 

For simplicity, we assume that the couplings are bidirec- 
tional, which means the corresponding graph is undirected, 



and K; 



K 



K.k with k denoting the k link 
We also assume that each coupling gain is symmetric 
positive definite, that is, = > . Assume now that 
the uncoupled dynamics f(xj,£) = f(xj,a, t) contains a 
parameter set a which is unknown to an arbitrary node <;. 
We then use the estimated parameter a in c;, with an adaptive 
law based on the local "coupling forces" 



a = PW T (x„t) K * (xj-x,) 



(5) 



where P and W(x cr ,t) are the same as those defined 
in ©. The network structure is illustrated in Figure [5] Note 
that the adaptation only uses feedback from the connected 
neighbors. 

To prove convergence, we define a Lyapunov function 



V = - ( x J L K x 



i T P 



i~ 



a) 



where Lx = D/CD T is the weighted Laplacian matrix, 
and each weight corresponds to the coupling strength of 



4 adaptation 



Fig. 3. Coupled network with adaptation 



that link. is symmetric positive semi-definite since K. is 
symmetric positive definite. We show that 



V 



x T L K x 



x T L K ( 



a T P *a 

f(xi,a, t) 

f(x„,a, t) 



- L^x ) 



= x T D/CAD T x - x T hi 



= x T ( L KA - L 



where 



L KA - D (tCA) s D 1 



The matrix AeK T 
fc th diagonal entry 

[A] fc = 



is a block diagonal matrix with the 



1 df 

o dx 

fch 



(xj +x(x. ( -Xj-)) rfx 



corresponding to the k link, which has been assigned with 
an arbitrary orientation by incidence matrix D, for instance 
from node i to node j. The matrix (JCA) S , the symmetric 
part of /CA , is also a block diagonal matrix with fc th 
diagonal entry [(/CA) s ] fe = (K fe [A] fe ) s . 

Lemma 1: The matrix L^a — L^; i s negative semi- 
definite if 

A^ +1 (L/c 



> max X m ax 
k 



(K fc [A], 



(6) 



,TlT 



For 



A„(L) 

See Appendix I VI-BI for the proof. 

Lemma 2: Given any vector x = [xf , x^, . . . , x^ 
a coupled network, if the condition (|6j is true, 

x T ( L KA - l£ ) x = 

if and only if xi = X2 = • • • = x„. See Appendix IVI-CI 
for the proof. 

Theorem 4: For a connected dynamic network (0 with 
adaptation l|5) adding to an arbitrary node <;, the states of 
all the elements will converge together asymptotically if the 
condition is true and all the states are bounded. 

Proof: Similar to the proof of Theorem|5] if all the states 
are bounded, we can conclude the boundedness of V, which 
then leads to the asymptotically convergence of V to zero 
if the condition (|6) is true. With Lemma |2] this implies 
immediately that xi = X2 = • • • = x n . □ 

Remarks: 

• Condition (|6j in Lemma [0 is in fact very similar to 
those given in Theorem|2] If X max (Kk[A}k) s is positive, to 



guarantee synchronization one needs a connected network, 
an upper bounded X max (Kk[A]k) s , and strong enough cou- 
pling strength JC. This result thus implies that adaptation 
will not significantly change the network's synchronization 
ability. 

Assuming m = 1 and all the coupling strengths = 
k , a sufficient synchronization condition for a coupled 
network without adaptation is 



while it is 



K 2 > 



A 2 (L) 
A n (L) maxiJi 



A 2 (L) A 2 (L) 

with adaptation. Here J JS = (df/dxi) 3 . 

• When the coupling gains are only positive semi- 
definite, extra restrictions have to be added to the uncoupled 
system dynamics to guarantee globally stable synchroniza- 
tion, similarly to the fixed-parameters result in [26]. See 
Appendix IVI-CI for details. 

• Theorem [4] requires the states to be bounded. Bound- 
edness can be shown following the same steps as we did in 
Section [Hi] In fact, since V < 0, we know that a and Vfc, 
Xft = Xj — Xj are bounded. Thus the boundedness of the 
states are determined by the Input-to-State Stability of the 
system 

Xj = f(Xj,t) +u, : 

Convergence of the estimated parameter set a can also be 
concluded with the same analysis as that in Section HTll 

• The result of Theorem |4] still holds if the parameter set 
a is unknown to multiple nodes and adaptations are added 
to these nodes simultaneously. This can be shown using the 
Lyapunov-like function 

1 

- ( x J L//c x + ^ ■ 



V = ^ (x T L K x + Vafp-% ) 



The nodes holding the real parameters can thus be con- 
sidered as "knowledge-based" leaders, which are different 
from usual "power-based" leaders [8], [11], [22], [26], 
which specify desired trajectories for the network by uni- 
directionally coupling to it. At the limit all nodes could 
be adaptive, although they may then converge to any odd 
parameter set — while all states will converge together, the 
desired individual behaviors (such as oscillations) may not 
be preserved depending on initial conditions. Note that both 
"power" and "knowledge" leaders may be virtual. 

• Note that in Lemma |2] the Xj may actually tend to 
zero together. We should exclude this possibility, and the 
possibility that any component of x, converges to zero, with 
dynamic analysis — for instance by showing that zero is an 
unstable state. 

Example 4.1: Consider a group of FN neurons connected in a 
general network 

(Vi = Vi(a — Vi)(vi — 1) — Wi + I + 22 kij( y j — y i) 
Wi — /3vi — -ywi i = 1, . . . , n 



Assume that all the parameters are now unknown to the node ?. 
We add adaptation 



where 



a = PW J 
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with W the same as that defined in Example 13.11 Simulation 
results are illustrated in Figure El □ 



Fig. 4. Simulation results of Example 14.11 The network contains four 
FN neurons connected in a two-way ring as Figure [3] The real parameters 
and the coupling gains are the same as those in Figure |2] The matrix 
P = diag{0.06, 3, 0.0002, 0.04} . All the initial conditions are chosen 
arbitrarily. The plots are (a), states Vi versus time, (b).V versus time where, 
V= \ (x T L K x + a^p-'a). 



B. Leader Combination 

We have mentioned that there can exist different leader 
roles in a network, ones with power and ones with knowl- 
edge. And in fact, these different types of leaders can 
co-exist. A leader guiding the direction may use state 
measurements from its neighbors to adapt its parameters 
to the values of the knowledge leaders. 

Consider such a leader-followers network 

x = f(x ,t) 

Xj = f(Xj,t) + ^ Kjj(xj - Xj) + 7;K ;(xo - Xj) 

where i = 1, .... n, x is the state of the group leader, 
7i = or 1, and Mi does not include the links with x . 
Adaptation can be added to any node(s) inside the network, 
such as 

n 

a = P W T (xo,t)^ 7i K oi (x ; -x ) 

i=l 

for the leader, or 

a ; = PiW T (x i7 t)(^ Kji (x j -x i ) + 7i K oi (xo-Xi) ) 

jeJVi 

for the followers. To prove convergence, we define several 
Laplacian matrices 

• Li/c, the weighted Laplacian of the followers network 
(excluding both the leader and the links from the leader). 

• Lx, the weighted Laplacian of the leader-followers net- 
work, which is not symmetric since we have uni-directional 
links. Moreover, 



-•K — 





-b C 



diag{7i K 0i } 



Note that C is symmetric positive definite if the whole 
leader-followers network is connected. 
• L/C, the weighted Laplacian of the leader-followers net- 
work if we consider it as an undirected graph. Thus, 



Lr — Lv 



b 



Define 



V = - ( x J L K x 



3 



a-P"^ ) 



where j may be equal to any number from to n. We have 

f(xi,a, t) 



V 



x L/c ( 

f(x„,a, t) 
= x T ( L KA - L£ L k ) x 



See Appendix I VI-EI for the conditions for L^a — L^J Ljc 
to be negative semi-definite. Following the same proofs as 
in previous sections, this then implies that all the states x^, 
j = 0,1, ... ,n converge together. Parameter convergence 
conditions are also the same. 

V. Concluding Remarks 

Coupled networks with adaptation are studied in this 
paper. We showed that synchronized behaviors can be 
achieved under similar conditions as we derived in [22], 
[26], which also guarantees parameter convergence if the 
stable system behaviors are sufficiently rich. Different kinds 
of leaders may coexist in the network. Current work in- 
cludes the applications of adaptive network, and the inves- 
tigation of coupled networks with distributed controllers. 

VI. Appendices 

A. Boundedness of Coupled FN Neurons 

For notation simplicity, define u = I + k(v2 — t>i) and 
w = w/ \/~j3. The dynamics of the first neuron changes to 



{vi = vi(a — vi)(vi — 1) — y/]3wi + u 
W\ — y/j3v\ — 7Wi 

Define U = | ( vj + u)f ), 

U = — (vi — a)(vi — \)v\ 



(7) 



UVi 



Since u is bounded, there must exist a large but bounded 
number vq > 0, V |vi| > v<j, U < 0. We denote the region 
\vi | < vq as Q. 



If the system {7) starts inside f2, since the dynamics of w\ 
is linear and strictly stable, vi and w\ are always bounded 
as long as the system stays inside f2. In fact, 



\Mt)\ 



< 



|(^(0) 
|^(0)| 



-ft 



,-7*) 



w (1 

7 



Thus, for any initial condition |u)i(0)| > vq^-, we have 
|wi (t) | < |iDi(0)|, which implies that the bound of 
inside f2 is max(«o— , |u>i(0)|). 

Suppose that at some moment, the system leaves fi 
through point (vq w ut)- Since U < outside f2, v\(t) 
and w\(t) will be both bounded by v$ + wl ut until the 
system trajectory re-enters Q, at which moment we should 
have | Wi„ | < I^Wtl- See Figure [5] for an illustration. The 
proof is similar if the system starts outside 51. 







1 * 


W 0Ut 






Q, 

























v o 


>- 



Fig. 5. Illustration of a solution trajectory of the system leaving and 
re-entering the region Q : \v±\ < vo. 

Thus, xi = [ Vi wi ] T is always bounded, which 
leads to asymptotic convergence of x to according to 
Theorem |3 Moreover, since the two FN neurons synchro- 
nize along a limit cycle, the convergence of W(x2,i)a to 
zero implies that of a. 

B. Proof of Lemma 1 

For notation simplicity, we first choose m = 1. Since 

Lka - tfc = D ( (/CA) S - /CD T D/C ) D T 

we know that is always one of its eigenvalues, with the 
corresponding eigenvector v = [1, 1, . . . , 1] T . Assume that 
the eigenvalues A^Lxa), Ai(Lje), ^ifiiKA — Lje) are all 
arranged in increasing order for i = 1, 2, . . . , n. According 
to Weyl's Theorem [7], for two Hermitian matrix A and B, 

A fc (A) + Ai(B) < Afe(A + B) < A fc (A) + A„(B) 

for each k ~ 1,2, ... ,n. Thus, we have 

A„-fe+i(L;cA — L^;) < A„(L/ca) — Afe(L^) 



which implies that, Vfc > 1, \n-k+i(f->K.A — Lj^) < if 



A n (L/cA) < A 2 (L^) 



Therefore, A„(Lxa — L 
negative semi-definite 

J K 



0, that is, L^ca — L£- is 



In fact, A 2 (L^) = A^Lyc). Assume 



maxA ma3; (K fe [A] fc ) s = A 

k 

If A < 0, we have A„(Lxa) < and both the conditions l|8} 
and are always true. If A > 0, then 

A n (L/c a) < A A n (L) 

and the result in Lemma ^ is concluded. 

In case rn > 1, we can follow the same proof except that 
zero eigenvalue here has m multiplicity, and the correspond- 
ing eigenvectors {vi, V2, . . . , v m } are linear combinations 
of the orthogonal set [I, I, ... , I] T where I e M mxm [ s 
identity matrix. 

C. Proof of Lemma 2 

For a real symmetric matrix, the state space has an 
orthogonal basis consisting of all eigenvectors. Without 
loss generality, we assume such an orthogonal eigenvec- 
tor set of L/ca - as {vi, v 2 , . . . , v m „}, where 
[vi, V2, . . . , v m ] = [I, I, . . . , I] T are zero eigenvectors. 
Therefore for any vector x, we have x = X^wi an d 

x T ( L KA - Lje ) x 



mn 

= E fcjV * T ( Lka ~ l k ) E fc * v 



i=l 



»=1 



mn 



kivj ( L KA - L| ) 



i=m-\-l 
mn 

E 

i—m-\-l 



mn 

E ■ 

i— m+1 



where A, is the i th eigenvalue corresponding to the eigen- 
vector Vi, and Xi < Vi > m for a coupled network if 
the condition l|6} is true. Thus, x T ( Lxa — L^- ) x = 
if and only if x = YnLi > tnat is> 

x = [1,1, . . . ,I] T [fci, . . . , k m ] T = [x ,x , . . . ,x ] T 

where x = [fci, . . . , k m ] T . 

D. Positive Semi-Definite Couplings 

Assume the coupling gain of the k th link is 



Kz 



Ki 




where Ki;. is symmetric positive definite and has a common 
dimension to all links. We divide the uncoupled dynamics 
J, and in turn the block diagonal entry of A into the form 



[A]* = s (*,t) 



Jfc(x,i) = 



Jn 
J21 



J12 
J22 



(8) 



where x is a value between the states of two neighboring 
nodes x, and Xj, and each component of Jk has the same 



dimension as that of the corresponding part in K^. Re- 
define the function V as 
1 



V 



( x t Lk;x + x T L Y x + a T P x a ) 



where Ly = DYD T is a weighted Laplacian based on 
the same graph as Ljc but different weights 


K 2 

Using a modified adaptive law 

a = PW T (x 5l t)^(K + Y) j5 (xj-xj 

we can show that 



V = 



( L (/c^ 



Y)A 



X L(£ +Y )(A-A) X 



where we define [A]fc = 

M fc = 

[(/C + Y)(A-A)] fc = 



Jn 





x T ( Wa - 
so that 



T 2 



KiJn 





Ki(Jii-J u ) 
K 2 J 21 



KiJia 
K2J22 



A non-positive V" can be guaranteed if 



• L 



< , which can be satisfied under similar 



condition as (jfji; 

• (£ + Y)(A-A) < 0, which is true if V k, the 
symmetric parts of Ki(Jn — Jn) and K2J22 are both 
negative definite, and <T max (Ki3i2 + J^Kj) is bounded. 
An explicit condition can be derived with feedback combi- 
nation analysis [22], [26]. 

The rest of the convergence proof are the same as that 
of positive definite couplings. 

E. Leader Combination 



Similarly to the proof of Lemma IVI-BI choose m = 1 
for notational simplicity. It can be shown that is always 
one of the eigenvalues of L^a — Ljc, which is negative 
semi-definite if 

A„+i(Lxa) < A2(L^ L/c) 

and its only eigendirection for the zero eigenvalue is v = 
[1,1,..., 1] T . Since 

b T b -b T C 
-Cb C 2 



we have 



A 2 (L£ L k )>A 1 (C 2 ) = A 2 (C) 



from the interlacing eigenvalues theorem for bordered ma- 
trices [7]. Thus a sufficient condition for synchronization 
is 

Ai(C) > A„ + i(L/ca) 



This condition is equivalent to the three requirements we 
listed in the first remark in Section I1V-AI Note that the 
connectedness condition refers to the whole network, while 
the subnetwork containing only the followers may not be 
connected. 

If all the coupling strengths are identical with gain k, the 
synchronization condition is 

A„ + i(L) maxjJis 



k > 



A 2 (L + diag{ 7 J) 



where L is the Laplacian matrix for the subnetwork contain- 
ing only followers, and L for the whole undirected leader- 
followers network. 

The proof is similar for m > 1. 
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